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Abstract 

> , 

CN ! A temporal analysis of the noise is performed, and non linearities 

are taken into account. We then extend the correlation method to 
■^j- ■ groups of several pixels to derive the interpixel capacitance of a de- 

tector, found to be x = -0.0263 ± 0.0020 (stat)±0.0040 (syst.) All 
measurements are consistent to a sub-percent accuracy. 



X ■ 1 The Apparatus 



The measurements described in this paper were carried out in a dedicated 
setup built to evaluate Hawaii detectors H2RG from Teledyne (ex Rockwell). 
The detector was on loan from LBNL in view of the evaluation of its perfor- 
mance when used in a spectrograph for the JDEM project pp. The cryostat 
can be operated in a range of temperature extending from 100K to 160K 
with fluctuations smaller than 0.1 K, and its equilibrium temperature in the 
absence of heating is about 110K. It was designed so as to ensure a variation 
rate of temperature smaller than 0.5 K/minute whatever the liquid Nitrogen 
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flow. A mirror located in the cryostat ensures a uniformity of illumination 
better than 1% over the full detector in the tests described in this paper. 

The readout cards were adapted from the acquisition system of the OPERA 
neutrino experiment. Their schematic layout is shown in Figure 1, and their 
configuration on the cryostat in Figure 2. This acquisition system was actu- 
ally used in the tests performed on the SNAP Spectrograph demonstrator in 
the Near infrared range [2J [3] . 



2 The calibration scheme 

We adopt the standard method of calibration as in [HE] taking advantage of 
the stochastic Poisson fluctuations from frame to frame under illumination 
by a Light Emitting Diode. The overall variance will be the sum of the con- 
tributions from readout noise, common mode noise, and stochastic noise, the 
latter being easy to extract since it is the only one to depend on illumination. 
The H2RG detectors have been extensively tested in the SNAP context as 
reported in [6], [10] . These authors also investigated the spatial correlations of 
the signal noise. The method proposed here is a variant where the emphasis 
is on a redundant determination of the interpixel capacitance using groups 
of pixels, and a temporal analysis of the LED signal is performed, as would 
occur in an actual flux measurement. Non linearities are taken into account. 

The ADC responds to a voltage change between two frames according to 
AA = kAV where the calibration coefficient k is the product of the emitter- 
follower factor ( 0.85) and of the ADC conversion factor (70 fiV/9.5/ADCU), 
and V is the pixel voltage. 

For each LED setting, 7 to 10 consecutive measurements were performed. 
The first two were ignored as the detector is not in a stationary regime, 
and the last five were averaged to obtain our results. The spread between 
consecutive exposures with the same LED intensity was used to estimate 
the measurement errors: their origin is NOT statistical, NOR due to a slow 
variation of the LED intensity. 

In a single pixel detector, the voltage change is linearly related (in the ap- 
propriate conditions) to the number of electron stored as AV P i X = qAN e /Co 
where AN e is the number of electrons stored in the pixel, q the electron 
charge, and C the capacitance of a single pixel, quoted by Teledyne as 40ff. 
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The conversion factor f e = e/ADCU is then expected to be 

AN e C 

fe = -TT = 7^ ~ 2 - 15 
AA kq 

The stochastic fluctuation of the charge on the pixel is SQ = q5N e , and 
the variance (time average) of the ADC readouts between two frames arising 
from Poisson fluctuations is < (5A) 2 >= (kq/C ) 2 < 5N 2 >= (kq/C ) 2 AN e 
Where AiV e is the accumulated number of electrons between the two frames 
considered (to be distinguished from the fluctuation SN e of this number). 
We can substitute AN e = AQ/q = C AV/q = C AA/(kq) which yields 

(SA) 2 = ^AA = l/f e AA (1) 

The conversion factor f e = e/ADCU is the inverse of the slope in the relation 
between the variance of the stochastic contribution to the noise (in ADC 
units) in a given time, and the flux accumulated during the same time (also 
in ADC units). The capacitance Co of the pixel can be obtained once k has 
been measured. 

In a multipixel detector, we adopt a description close to the one proposed 
by |12] , but more general. The charge in pixel i is now also dependent on 
the voltage in pixel j, and the relation is given by the electrostatic influence 
matrix. AQi = CijVj (summation is implied). The ADC response in pixel i, 
AAi will be given by: 

AAi = kqUl^A.X, 

where ANj is the change in the number of electrons in pixel j. The stochastic 
noise of pixel % can again be derived from the Poisson fluctuations in the pixels 

(6 t ) 2 = {kqfiC^iC^HN^N, 

If diffusion is negligible, there is no correlation in the numbers iVj for different 
pixels, and the time average of the product is < SNjSNi >= SjiANj (where 
Sji is the Kronecker matrix, and ANj the nb of electrons collected in pixel j 
between consecutive frames. Using ANj = l/(kq)Cj m AA m we then find 

(5 l ) 2 = fcg^(C) ij . 1 C^Q m AA m 

j 

So that 

(Si) 2 = kq J £C^ 2 C jm AA m (2) 
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This general expression will now be simplified by the use of convenient ap- 
proximations of the matrices C and C _1 . If we consider only the coupling of 
adjacent pixels, with x = Cj/Co (x is negative) in a homogeneous detector, 
as in [3], the influence matrix C for a group of 5 pixels centered on pixel 
i = as shown in Figure 3b) is C = CqM, where the (5 x 5) matrix M is of 
the form 



M 



This matrix can be inverted exactly, but as the electrostatic coupling ratio x 
will be found to be small, the inverse matrix can be obtained by substituting 
— x to x when the second order terms x 2 are neglected. 

We have now obtained for the single pixel noise of a homogeneous detector 
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kqC^SNj 



k-^r(SN - arfJVi - x5N 2 - x5N 3 - x5N A ) 



evaluating the variance from both sides of this formula, we find: 

< W 



> 



^) 2 (AiV )(l + 4a; 2 ) 



a 



As the illumination is uniform, all fluxes are equal and the mean flux in pixel 



is AN 



<(5 



Co 



AAr 



> 



kq 



C (l-4x) 



feAA 

AA = siAA 



(3) 



Under a uniform illumination, the conversion factor (e/ADCU) was however 
found to be 

, AiV C 

Je 



AA kq(l - Ax) 



so that 

fe 



Si(l - Axf 



(4) 



As x is negative, the conversion factor will now be smaller than the inverse 
of the slope, and smaller than expected for a 'single pixel'. The edge effect 
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from adjacent pixels is expected to decrease for larger groups of pixels, and 
the relation between slope and conversion factor should be closer to the single 
pixel case. 



3 Noise fluctuations for groups of pixels 

The previous formulae allow us to measure the ratio x = Ci/Cq by comparing 
the relation between noise and flux for groups of pixels, as this changes the 
weight of the contribution of the neighbouring pixels. 

2 pixel groups (as shown in Figure 3a) 

The 8x8 capacitance matrix C to be considered follows from Fig. 3 a): 
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For a typical value of x = —0.02, the diagonal elements of the inverse ma- 
trix vary from 1.0004 to 1.0016, while the off-diagonal elements of adjacent 
pixels differ from 0.02 by less than 510 -5 , the other coefficients are smaller 
than 410~ 4 . Neglecting all small offsets, we find: 



^(SQ^l - x) + 5Q 2 (1 - x) 



) 2 = (^n2-4x + 8x>)~A^ 



and substituting for the expression of A7V in terms of the observed ADC 
shift AA in 1 pixel between frames 
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AN <> = 5<fe) AA » 

5 pixel groups (Figure 3b) 

We find in the same way: 

k= ^(5Qo(l-4x) + (l-x)E i=M =4^ 

— 2x Y^i=5,i=8 $Qi ~ x Si=9,i=12 &Qi) 

So that 

(5 5 ) 2 = (^) 2 (5 - 16a; + 40a; 2 ) ~ (^) 2 AiV (5 - 16a;) 

and using as before the relation between ANq and AAq, the ADC shift: 
,„ 2 / kq . 5 - 16a; + 40a; 2 . 

9 pixel groups (Figure 3c) 

and for 9 pixels forming a (3 x 3) square (see Figure 3): 
5 9 = ^Wo(l - Ax) + (1 - 3x) E t =i,7 H d) «?i 

+ (1 — 2x) X)i=2,8(euen) $Qi ~~ x 1^9,20 ^Qi) 

(£ 9 )2 = (M^) 2 (9 - 48x + 80x 2 )AiV 
Co 

r9 fco 9 - 48a; + 80a; 2 . , 



4 The measurement method 

The cryostat temperature was set at 110K, the equilibrium value in the ab- 
sence of any heating, during these tests. In contrast with previous publica- 
tions, which used the spatial correlations under different illumination condi- 
tions ([I]), we have used the time variation of the signal to evaluate the noise 
and the correlation, as a training for the flux measurements anticipated in 
the near future. The frames are grouped into 'bursts', the number of frames 
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in each burst decreasing from 60 (LED current of 10 fiA), to 15 (LED cur- 
rent of 100 yuA). The variance of the differences of readouts for consecutive 
frames is then evaluated for each pixel, and each burst, and the measure- 
ment is repeated 'up the ramp' until saturation is reached. The results are 
then averaged over all pixels for each burst. At any LED setting at least 10 
exposures are taken to allow control of the error estimates. 

4.1 Non-linearities 

The distribution of the flux in ADC units averaged over all bursts for all pixels 
at a typical setting with a LED current of 40 /iA is shown in Figure 4 a). It 
is seen that in the absence of flatfielding corrections, its spread is 3.5%. This 
good homogeneity of the detector in the analysis window allows averageing 
of the measured properties over all pixels (the flux and variance of each pixel 
is still evaluated independently). A linear variation of dADC/dt and of the 
single pixel noise along the ramp is then observed in Figure 4 b). As the 
flux measurement is perfectly reproduced from one exposure to the next, the 
flux variation is an indication of the non-linear response of the system. The 
most likely source of this behaviour is the output FET, and the observations 
suggest a decrease of the transconductance as the grid voltage decreases with 
an increasing number of trapped electrons in the pixel well. The slopes of the 
variance and of the flux (normalised to their value at an output ADC value of 
19000 ADCU) are respectively a var = 0.000388 and a flux = -0.000336, and 
they correspond to similar (and opposite) relative variations, as expected if 
the equivalent thermal resistor is the inverse of the transconductance. The 
effective conversion factor given by the ratio variance/flux cannot be obtained 
without further corrections. All the following analyses have been performed 
with 2 different extrapolations, namely to the middle (32500) and to the lower 
values (10000) of the ADC dynamical range. To crosscheck the validity of the 
corrections, different LED illuminations will be compared. The conversion 
factors found are expected to differ by 5% as a consequence of the non linear 
response, but the interpixel capacitance obtained should be the same up to 
systematic errors. 

4.2 The 2 pixel pairs 

It is intructive to consider separately horizontal (readout direction) and ver- 
tical pairs of pixels. Horizontal pairs are read consecutively, while vertical 



7 



Table 1: Readout Noise Horiz./Vert. pairs 



Pixel 


var(ADCU) 


group 




1 


121. 02± 0.77 


2-H 


249.17± 1.50 


2-V 


261. 24± 1.79 



pairs are separated by the time interval needed to complete the readout of 
the intermediate pixels. The variance quoted in Table 1 is obtained in the 
absence of any illumination by summing 2 pixels, measuring the difference 
between consecutive frames, estimating its variance for each burst (of typi- 
cally 60 frames), averaging over all bursts, and then over all pairs. The result 
for the mean variance is given. 

It is seen that intrinsic voltage changes in the chip during the readout 
impact the variance a level of about 4%. Analyse is in progress to shown 
soon results including a temporal common mode correction by subtracting a 
reference channel. 



4.3 The 5 and 9 pixel configurations 

The measured values of the variance and the flux for each pixel grouping is 
shown in figure 5 for 12 illumination conditions. The expected linear relation 
(after correcting for the non-linear response!) is observed in Figure 5 for the 
data and the slopes found from the extrapolation to lower flux values are 
given in Table 2. 

The ratio x of the interpixel capacitance to the pixel capacitance Co is 
found from the comparison of the slopes for 1,2, 5, and 9 pixels: To modify 
the impact of the fluctuations from adjacent pixels, we now evaluate their 
contribution to the variance of larger groups of pixels. The interpixel capaci- 
tance x = C c / C can then be derived for each pixel grouping from equations 
(5),(6),(7). 
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Table 2: slopes 



pixel 


Slope 


group 




1 


0.4021 ± 0.0026 


2 


0.8458 ± 0.0049 


5 


2.2139 ± 0.0218 


9 


4.2095 ± 0.0550 



4.4 Conversion factor and effective Interpixel capaci- 
tance 

The values of the interpixel capacitance found in all cases are given in Table 
3 with the linearity correction extrapolated to the lower ADC range (10000). 
They are compatible, with small residual systematic shifts, and we shall 
perform a weighted average. 



Table 3: The Slope ratios and Interpixel capacitance 



Group 


ratio 


X 


s2/sl 


2.1028 ± 0.0183 


-0.02455 ± 0.00416 


s5/sl 


5.5063 ± 0.0651 


-0.02947 ± 0.00352 


s9/sl 


10.4696 ± 0.1529 


-0.02916 ± 0.0029 




x = 


-0.0282 ±0.0020 



The value x = —0.0282 ±0.0020. The same averages performed at ADCU 
= 30000 leads to x = -0.02440 ± 0.001 For our final result, we take the 
average of the two estimates and attribute a systematic error of 0.0040 to 
account for the small difference in the non-linearity corrections (mid range 
and low range). 

C 

x = mt = _ 02Q3 ± o.0020(stat) ± 0.0040(syst) 
Co 

The (single pixel) conversion factor in the lower ADC range can now be 
obtained as 

AQo 1 
h ' AADCo Sl (l-Ax) 2 
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2487 2 ' 487 2. 03H e/ADCU) 



(l-4x) 2 1.221 

and it is 10% larger at full well. The conversion factors found in the 2 ex- 
trapolation methods differ by 5%, as expected from the nonlinear behaviour 
seen in figure 4b): they are evaluated at different values of the ADC range. 
The result is very close indeed to the value of 2.15 derived from an assumed 
pixel capacitance of 40ff in section 2. 



4.5 Diffusion and actual Interpixel capacitance 

The effective value of the interpixel capacitance found would NO The the true 
value if 'fast' diffusion would occur as suggested by [5J: while the interpixel 
capacitance increases the fluctuations, diffusion from one pixel to the adjacent 
ones would reduce them. A detailed computation shows that the contribution 
to Xd from diffusion is equal to the fraction fd of electrons migrating to the 
adjacent pixel for 2 and 5 pixel groups, but is only Xd = /d/3 for the 9 pixel 
group. Given the consistency of the previous results between 9 and 5 pixel 
groups, we obtain f d < 0.037 (3 a limit). 



5 Conclusions 

A general correlation method has been proposed, which allows strong cross- 
checks for internal consistency of the observations between different pixel 
groups. We have shown that the non-linearities which are seen are consis- 
tent with the effect of a transconductance variation in the output FET of 
the detector, and that they can be corrected to a sub-percent accuracy. The 
remaining systematic errors are in the 10 _3 range. It has also been shown 
that the use of the reference channel have to studied in that frame of work. 
The imapct of that studies will be shown in a forthcoming paper. The cal- 
ibration of the readout set-up will be studied in order to describe the noise 
performance achieved under different conditions. 
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Figure 1: Scheme of the acquisition card 
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Figure 2: Cryostat with the readout cards 
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gure 3: Map of the 2,5, and 9 pixels groups 
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Figure 4: a) flux distribution ADCU/pixel/frame for an LED setting at 40 
fiA b) variance (ADCU) 2 and flux (ADCU) as a function of the ADC value, 
for a LED setting at 20 fiA. \q 
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Figure 5: noise(without reference subtraction)/Flux for single pixels and 
groups of 2,5 and 9 pixels 
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